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Abstract 

We review curvature-based hyperbolic forms of the evolution part of the Cauchy 
problem of General Relativity that we have obtained recently. We emphasize first 
order symmetrizable hyperbolic systems possessing only physical characteristics. 
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1 Introduction 



Application of general relativity to fully three dimensional problems in astrophysics and 
cosmology has provided a driving force motivating further study and reformulation of 
Einstein's equations in 3+1 form. The vital role of the four constraint equations in 
setting up initial data for numerical evolution and in carrying out mathematical studies 
of gravitational field configurations has long been recognized. It led in the 1970's to a 
general, mathematically rigorous, and useful formulation of the constraints. (See, for 
example, the review The standard evolution equations for the spatial metric g and 
the extrinsic curvature K were already known from a straightforward decomposition of 
the spacetime Ricci tensor (cf. the cases of zero shift ||^, arbitrary shift 0, and synthesis 
from the spacetime viewpoint and further developments |^). It was widely believed that 
these equations (first order in time derivatives; second order in space derivatives) were 
adequate for applications. However, they do not constitute a hyperbolic system leading 
to a proof of causal evolution in local Sobolev spaces of g and K into Einsteinian 
spacetime. Hence, further study of evolution equations is essential. 

In this paper, we review methods that we have used recently to obtain hyperbolic 
systems for the evolution of (g, K) with only physical characteristics and which prop- 
agate curvature along the physical light cone. We emphasize first order symmetrizable 
hyperbolic systems. 

2 Einstein-Ricci Hyperbolic System 

What we may call the "Einstein-Ricci" (ER) system is a spatially covariant hyperbolic 
formulation of the Einstein equations constructed from first derivatives of the spacetime 
Ricci tensor. The ER equations are obtained from the 3+1 form of the metric 

ds'^ = -N^dt^ + gij{dx' + (3'dt){dx^ + ^^dt), (1) 

where > is the lapse, /3* is the shift, and the properly Riemannian metric gij is 
the spatial metric of a spacelike slice S. (S is understood to be a generic spacelike 
leaf of a foliation of a globally hyperbolic spacetime: locally, t=const.) It is convenient 
to work in the cobasis 6^ = dt, 6^ = dx' + (3'dt, with = A 6^ and only 

C'oj = —C^jo = djP^ 7^ 0. The connection one-forms uj°-p = "^"p-yO^ are given by the 
connection coefficients 
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where T denotes a Christoffel symbol. In particular, the extrinsic curvature of S is 
Kij = —Nj^ij, or, upon evaluation of 7°jj, 

dog^j = -2NKi,, (3) 

where do = dt — £13 evolves t-dependent spatial objects in the direction orthogonal to 
S, and £/3 is the Lie derivative in E along the shift vector. The relation (^) serves as 
the evolution equation for gij. Letting an overbar signify spatial objects, we note for 
future reference that (jij = g^j, g^^ = g^\ and Yjk = r*jfc = ^^jk- For completeness, we 
also state 7^00 = NV'N, 7%o = N-^doN, 7%, = 70^0 = N-^ViN, 7^, = -NK^, and 
7^io = I'^oi + di(3^. 

The 3+1 decomposition of the spacetime Riemann tensor is given, for example, in 
Ref. [4]. Likewise, the Ricci tensor is given identically by 

Ri, = Rij-N~^doKij+HKij-2KikK''j-N-'ViVjN, (4) 
i?o. = NV'{Hg,,-K,^), (5) 
Roo = NdoH-N^K^kK"'' + NWkN, (6) 

where H = K^k = K. The heart of the ER system is defined by 

Vtij = doRij — 2V(j-Rj)o. (7) 

Upon working out this identity using (|)-(||) and substituting the Einstein equations 
Ra/3 = Pa/3 = 87r(TQ,^ — ^ga/sT'^-y), with G = c = 1, wc obtain the equation 

- NnKij = Jij - Sij + (8) 

where the physical wave operator for arbitrary shift is □ = —N~^dQN~^do + V^Vfc. We 
study the vacuum case flij = here. In general, the value of Qij is defined by matter 
fields. 

The detailed form of the right hand side of (|]) can be found in Refs. [5]- [7]; the present 
conventions are those in Refs. [6], [7]. Here we point out that Jij = Jij[2 in g; 2 in A^;l in K]. 
(The numbers in the bracket indicate the highest order derivatives that occur.) 

The slicing-dependent term S^j is given by 

Sij = -N-^ViVjidoN + N^H). (9) 
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We observe that the term ViVjH = g"^^Vi'VjKmk would spoil the hyperbolicity of (|); 
therefore, Sij must be set equal to a functional involving fewer than two derivatives of 



Kij. Notice that 



^oA^ + N^H = -iV='7° = (7^/2^(^-1/2 ^ g^^doa, (10) 



where g = det gij > 0, 7° = g°'^'j'^ai3, and a = g^^/'^N is a lapse function of weight 
-1, first introduced in the "algebraic gauge" and subsequently as a "proper-time 
gauge" IP and in connection with the "new variables" program for general relativity 
in Hamiltonian form[]10|, where a = iV is called the "densitized lapse." The algebraic 



gauge has a simple relationship to "harmonic time slicing," 7*^ = 0, as we see in ([T0|). 

One most easily deals with Sij by freely specifying a(t, x) > 0. (The shift (3\t,x) 
is also arbitrary.) A "gauge source" f(t,x) can be used, as observed by Friedrich (see 
TT[| for references); here, f{t,x) = c^ologa. Then the first term in ([lOD yields a general 



harmonic time slicing equation for A^, namely 

doN + N^H = Nf{t,x), (11) 

or equivalently, 

doig~'^'N)=g-'/'Nf{t,x). (12) 

In this case, the system composed of (|l^) and the third-order equation resulting from 
combining (^ and (|^), whose principal term is ndogij, is quasi-diagonal hyperbolic |]12 



for the unknowns g and g~^^'^N. Therefore, the "third-order" ER equations have a 
unique, well-posed solution in a suitable Sobolev function space (local in time and in 
space; global in space can be obtained as well), given appropriate initial data |13|. 

We wish to verify that every solution of the third order ER system, equivalently {(§), 
(I), (|12D or (|ll|)}, is a solution of the Einstein equations given suitable initial data. Thus, 
suppose g, K, and satisfy these equations. Then, by using the contracted Bianchi 
identities V^G/3q, = 0, and these ER equations, it follows that G/^a = if g and K 
satisfy initially the usual momentum and Hamiltonian constraints (Cj = —2N~^Roi = 0, 
C = 2G^o = 0), if doKij satisfies Rij = initially, and if initial data {a > 0,/3*} 
or {N > 0,(3^} are given |^, ^. Conversely, because every globally hyperbolic metric 
solution of Einstein's equations can be given in a harmonic time slicing |^], it follows 
that all globally hyperbolic solutions of Einstein's equations can be reached by solving 
this Einstein-Ricci system. 

The ER system as given holds for any spacetime dimension. For four dimensions, 
it can be written in first-order symmetrizable hyperbolic (FOSH) form. To put a sys- 
tem in first-order form, one introduces auxiliary variables in place of derivatives of the 



4 



fundamental variables and adds additional equations to describe the evolution of these 
new variables. The crux is whether this process terminates. Generally existence of a 
higher order wave equation on the fundamental variables is necessary to halt the pro- 
cess. Though we have a wave equation for K, presence of second spatial derivatives of 
g and N in Jij prevent the reduction to first order form from being obvious. However, 
the second derivatives of g occur only in Rij and Rijki terms; the latter, in three space 
dimensions, can be reduced to Rij, which in turn can be eliminated from Jij by using 
(I) . First derivatives of g are handled by an evolution equation [|], ^ for the Christoffel 
symbol f (inadvertently omitted in [6]). 

On the other hand, to handle the second derivatives of N, one first finds a wave 
equation for N. Then, by applying Vj to it, a wave equation for ViN, or Oj = Vj log A^, 
is obtained. From this point, reduction to first order form is straightforward [The 
wave equation for follows from applying Bq to (|TTD to get an equation with dodoN and 
doH. Using (^), we then find a nonlinear wave equation that can be written in terms of 
□A^. Then ViUN gives the wave equation naj. See Refs. [5], [7] for details.] One then 
shows that the first order system is symmetrizable and hyperbolic. 

In this system, the spacetime metric {N, Qij) evolves at zero speed along the direction 
do orthogonal to t=constant. Quantities with dimensions of curvature propagate at speed 
1 ("c") along the physical light cone. There are no unphysical speeds or directions as 
there are in many of the FOSH systems that have so far appeared in the literature |Tl| , 
15|. Of course, we except the FOSH form of the ER system (described here) and 



the mathematically equivalent, but more transparent, "Einstein-Bianchi" [T^, (EB) 
system presented in the next section. Having only physical characteristics is essential for 
an ideal match of physics and mathematics. In practice, it is also useful: for example, 
consider obtaining the gravitational radiation content of a numerical Cauchy evolution 
on a finite grid. The radiation is propagating at light speed and is either "extracted' 



or "evolved to null infinity" |TP| based on field values at a finite distance from the source. 
We know that gravitational radiation is curvature propagating at light speed and that 
is what the ER system describes. The description is even more explicit and transparent 
in the EB system. 



A "fourth order" ER system also exists |^0| . One forms the expression 



doQij + ViVjRoo (13) 

to obtain it. It has a well-posed Cauchy problem and is hyperbolic ("non-strict" [^) 
for any A^ > and It has been applied to the non-linear perturbative regime of high 
frequency wave propagation ^ and has been used to obtain an elegant derivation 



of gauge- invariant perturbation theory for the Schwarzschild metric [^| 
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3 Einstein-Bianchi Hyperbolic System 



The "Einstein-Bianchi" (EB) system discussed next was given first in [16] and with 
energy estimates and full mathematical detail in [17]. It is similar to an analogous 



system, obtained by H. Friedrich [|TT], that is based on the Weyl tensor and is causal 
but with additional unphysical characteristics. Recall that the Riemann tensor satisfies 
the Bianchi identities 

VaRf3-y,Xfi + ^'yRal3,Xfj. + ^ l3R-ya,Xfi = 0, (14) 

where we use a comma to stress the two separate antisymmetric index pairs {not to 
indicate partial differentiation). These identities imply by contraction and use of the 
symmetries of the Riemann tensor 

'^aR'^PM + Vm^/3A - VAi?/3^ = 0, (15) 

where the Ricci tensor is defined by 

R"'(3,afi = Rpfi- 

If the Ricci tensor satisfies the Einstein equations 

Rap = Paf3, (16) 

then we have 

VaR^f^M = ^Xpf^t^ - V;.P/3A. (17) 

Equation ([I^) with {aP'j} = {ijk} and (|17D with P = contain only derivatives of 



the Riemann tensor tangent to S; hence, we consider these equations as constraints 
("Bianchi constraints"). We shall consider the remaining equations in ( p!4D and (|1^) as 
applying to a double two- form A^/s^Xfj., which is simply a spacetime tensor antisymmetric 
in its first and last pairs of indices. We do not suppose a priori a symmetry between 
the two pairs of antisymmetric indices. These "Bianchi equations" are 

^oAhk,x^l + Vfcv4oh,AAt + V/iv4fco,A^t = 0, (18) 

VoA°,,A^ + V;,A\,Am = ^XP^^. " ^ ^^P^x = Ji^, (19) 

where the pair [X/i] is either [Oj] or [jl]. We next introduce following Bel p5|] two 
"electric" and two "magnetic" space tensors associated with the double two-form A, 
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in analogy to the electric and magnetic vectors associated with the electromagnetic 
two-form F. That is, we define the "electric" tensors by 




i,Oj - 



■0i,0j5 



(20) 



VihkVjlmA 



hk,lm 



while the "magnetic" tensors are given by 



(21) 



In these formulae, rjijk is the volume form of the space metric g. We note that: (1) If the 
double two-form A is symmetric with respect to its two pairs of antisymmetric indices, 
then Eij = Eji, Dij = Dji, and Hij = Bji. (2) If A is a symmetric double two- form such 
that Ao^p = A^a,x/3 = cgai3, then Hij = Hji = Bji = Bij and = Dij. Property (1) 
is obvious and (2) follows from the "Lanczos identity." ||2^ (We note that rjoijk = Nrjijk 
relates the spacetime and space volume forms in using the Lanczos identity.) 

In order to extend the treatment to the non-vacuum case and to avoid introduc- 
ing unphysical characteristics in the solution of the Bianchi equations, we will keep as 
independent unknowns the four tensors E, D, B, and H, which will not be regarded 
necessarily as symmetric. The symmetries will be imposed eventually on the initial data 
and shown to be conserved by evolution. 

We now express the Bianchi equations in terms of the time dependent space tensors 
E, H, D, and B. We also express spacetime covariant derivatives V of the spacetime 
tensor A in terms of space covariant derivatives V and time derivatives do of E, H, D, 
and B by using the connection coefficients in 3-1-1 form as given earlier. 

The first Bianchi equation (|1^) with [Xfi] = [Oj] has the form 



0. 



(22) 



A calculation yields (|2^) in the form 



do{r]\kHij) + 2NV[hEk]j + {Li)hk,j 



(23) 



(Li)mj = NK',r^\kHa + 2{y^hN)E^,+2Nr,\,K\kBt,^, 



(24) 
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The second Bianchi equation (|T^), with [A/i] = [Oj], has the form 



(25) 



where J is zero in vacuum. A calculation similar to the one above yields for the second 
Bianchi equation 

doE,^ - Nri^\\/hHij + (L2)i, = J.,o„ (26) 



(27) 



The non-principal terms in the first two Bianchi equations (23) and (|26|) are linear 
in E, D, B, and H, as well as in the other geometrical elements A^K and VA^. The 
characteristic matrix of the principal terms is symmetrizable. The unknowns -E'j(j) and 
ifj(j), with fixed j and i = 1, 2, 3 appear only in the equations with given j. The other 
unknowns appear in non-principal terms. The characteristic matrix is composed of three 
blocks around the diagonal, each corresponding to one given j. 

The j^^ block of the characteristic matrix in an orthonormal frame for the space met- 
ric g, with unknowns listed horizontally and equations listed vertically, {j is suppressed) 
is given by 







E, 


E2 


E3 




H2 


H3 




(26)i 




eo 















\ 









eo 










AT^i 


















-mi 







0)23 











eo 










(Hal 














eo 







(El 12 


I 


-iV6 


mi 











^0 


/ 



(2^ 



This matrix is symmetric and its determinant is the characteristic polynomial of the E, 
H system. It is given by 

-N'i^oe)iUn'- (29) 

The characteristic matrix is symmetric in an orthonormal space frame and the timelike 
direction defined by do has a coefficient matrix Tq that is positive definite (here Tq is 
the unit matrix). Therefore, the first order system is symmetrizable hyperbolic. We will 
not compute the symmetrized form explicitly here. 
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The second pair of Bianchi equations is obtained from (|T8|) and ( |T9| ) with [X^] = [Im]. 
We obtain from (|l^) 

doiv'hkV^ imDij) + 2Nr]^imV[kBh]j + {L^)hkM = 0, (30) 

{L3)hkM = 2Nr]'',^^K^i^r]\kDin + 2v'iU^ikN)Bh]j (31) 
+2NKi[hEk]„, + 2NKm[kEh]i + 2Hi^{y m]N)r]\k. 

Analogously, from (|I9]) we obtain 

Wln.B^J) - Nr]''\r]^imVhDkJ + {U)iM = -NJi^im, (32) 

iL,),,im = -N{TTK)r]^i^B,, + 2N7]'' ,^^Kh]Bih + 2NK\r]^ iraB^j (33) 

Consider the system (^) and (^) with [Im] fixed. Then j in riji^ is also fixed. The 
characteristic matrix for the [Im] equations, with unknowns and Bij, j fixed, with 
an orthonormal space frame, is the same as the matrix (^8]) . 

If the spacetime metric is considered as given, as well as the sources, the Bianchi 
equations (p3D , (p6|), (|30|) , and ( ^21) form a linear symmetric hyperbolic system with 
domain of dependence determined by the light cone of the spacetime metric. The coeffi- 
cients of the terms of order zero are VA^ or A^K. The system is homogeneous in vacuum 
(zero sources). 



4 Determination of (F, K) from Knowledge of the 
Bianchi Fields 

We next link the metric and connection to our Bianchi fields. This link uses and extends 
an idea introduced by Friedrich in his Weyl-tensor construction mentioned above. 



We will need the 3 + 1 decomposition of the Riemann tensor, which is [Q 

Rij,ki = Rij,ki + 2Ki[kKijj, (34) 

Ro^,Jk = 2NV[jKk],, (35) 

Ro^,OJ = N{doK,, + NK,kK^j + V,d,N). (36) 
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From these formulae one obtains those for the Ricci curvature given in Sec. 2: (|^), (^), 
and where, in this section, we will not denote K^j = TrK by H. The identity 
and the expression for the spatial Christoffel symbols give 

aof = V\NK,,) - 2V(,(iVir,/). (37) 

Therefore, from the identity (^51), we obtain the identity 

doT\, + NV^Ki, = K,,d^N - 2K\,d,)N - 2i?o(.,,/. (38) 

On the other hand, the identities (pB]) and (|^) imply the identity 

doKij + NRi, + VjdiN = -2NR\^Q, - iV(TrK)/Ci, + NRi,. (39) 

We obtain equations relating F and K to a double two-form A and matter sources 
by replacing, in the identities (||) and (|9D, -Ro{j,j)^ by (Ao{i,j)'' + ^''(j,i)o)/2, -R°i,oi by 
{APi^Qj + AP j Qi) /2, and the Ricci tensor of spacetime by a given tensor p, zero in vacuum. 
The terms involving A are then replaced by Bianchi fields E, H, D and B. 
The first set of identities (p8| ) leads to equations with principal terms 

doT\, + Ng^^d,,K,,. (40) 

To deduce from the second identity (^) equations which will form together with the 
previous ones a symmetric hyperbolic system, we set, we use algebraic gauge as in 
Sec. 2, by setting 

N = g^/^a (41) 

where a is a given positive scalar density of weight minus one, a function of {t, a;*). (Note 
that the present a is the of |jT6| and [0.) The lapse is now a derived quantity 
depending on g^^'^ and on a. The use of a, if T denotes the Christoffel symbols of g, 
imphes that 

T^h = d^hgN-di\oga. (42) 
The second set of identities ( pUf ) now yields the following equations, where N denotes 

doK,,+NdhT^, = iV[f'",,fV-(f^^ + '9aoga)(^'=,fc + c',loga)] (43) 
-N{did, log a - T^jdk log«) - N{Eij + Eji) - iV(TrK)Ki,- + Npi,. 
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The first set ( pS]) yields 



(44) 



We see from the principal parts of ( ^31) and (^) that the system obtained for K and 
r has a characteristic matrix composed of six blocks around the diagonal, each block 
a four-by-four matrix that is symmetrizable hyperbolic with characteristic polynomial 
-^^(^o^°)(^a^")- The characteristic matrix in a spatial orthonormal frame has blocks of 
the form 



Nil 



Nil 

io 





Ni2 



io 




^3 \ 




io J 



(45) 



5 Symmetric Hyperbolic System for (E, H, D, B, g, K, T] 

We denote by S the system composed of the equations (^), (p6D, (pOD, (|32|), (|^), (^), 
and (0), where the lapse function is replaced by g^^'^a. This system is satisfied 
by solutions of the Einstein equations whose shift (3, hidden in the operator do, has 
the given arbitrary values and whose lapse has the form N = (Clearly, any 

> can be written in this form.) From the results of the previous sections, we see 
that for arbitrary a and /3, and given matter sources p, the system S is a first order 
symmetrizable hyperbolic system for the unknowns (E, H, D, B, g, K, F). Note that the 
various elements E, H, D, B, g, K, and f are considered as independent. For example, 
a priori, we neither know that F denotes the Christoffel symbols of g nor that E, H, D, 
and B are identified with components of the Riemann tensor of spacetime. 

We now consider the vacuum case. The original Cauchy data for the Einstein equa- 
tions are, with (f) a properly Riemannian metric and ^|) a second rank tensor on an initial 
spacelike slice Sq, 

g|o = 0, K|o = V- (46) 
The tensors and i/' must satisfy the constraints, which read in vacuum, 

i?o. = 0, (47) 
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— Goo = -Roo — 2^ooR (48) 
= i?oo + 

with R = g^^Rij. The initial data given by <p determine the Cauchy data T^ij\o and thus 
Rij,ki\o- Then, Rij^ki\o is determined by using also i/'- To determine the initial values 
of the other components of the unknowns of the system S, we use the arbitrarily given 
data a and f3. In particular, we use = g^^'^a to find Roijk\o and to compute VjdiN\o 
appearing in the identity (|). We deduce from (^) doKij\o when Rij = 0, which enables 
Roi,oj to be found from (^). All of the components of the Riemann tensor of spacetime 
are then known on Sq. We identify them with the corresponding components of the 
double two-form A on Sq: the latter have thus initially the same symmetries as the 
Riemann tensor. We find the initial values of (E, H, D, B) according to their definitions 
in terms of A. 

Detailed existence theorems and proofs that solutions of the EB system, with initial 
data given as in the preceding paragraph are found in [18]. Here we give a less detailed 
argument, which is nevertheless completely rigorous. 

Let the initial data for the vacuum EB system be given as above. We know that our 
symmetrizable hyperbolic system has a unique solution. Because a solution of Einstein's 
equations in algebraic gauge [i.e., with given a{t,x)] — previously proven to exist |^ — 
together with its connection and Riemann tensor, satisfies the present EB system while 
taking the same initial values, that solution must coincide with the solution of the 
present EB system in their common domain of existence. 

The POSH form of the ER system and the EB system are completely equivalent 
mathematical systems, with the EB system perhaps having its unknowns arranged in a 
more transparent way. For instance, the characteristic Bianchi fields with respect to a 
fixed space direction that propagate along the light cone are the directionally transverse 
components of the Riemann tensor. For a detailed discussion of charateristic fields of 
the FOSH form of the ER system-also curvatures-see [7]. These fields for the EB system 
will be discussed in more detail elsewhere, as will "third order" and "fourth order" forms 
of the EB system, analogous to those of the ER system. 
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6 A Concluding Remark About Constraint Viola- 
tions 



It is well known that the standard 3+1 equations (0)-®, with time derivatives doQij and 
doKij only, were reformulated in terms of closely related variables and written in Hamil- 
tonian form in far-reaching work by Arnowitt, Deser, and Misner (ADM) and by 
Dirac ^3. Indeed, the standard 3+1 equations are often referred to as the "ADM equa- 



tions." However, it can be shown ||2^, |30[ that when there are violations of the constraints 
(as there generically are in numerical work), the standard 3+1 evolution equations and 
the ADM evolution equations are not equivalent. From an analytic standpoint, the 
standard 3+1 equations in terms of (g, K, Rij, Roi, G^o) produce a first order symmetric 
hyperbolic system for the constraint violations, to which suitable "energy" bounds for 
the growth of these violations can be applied, while the actual ADM equations in terms 
of (g, TT, Gfj_^) do not produce a hyperbolic system for evolution of constraint violations. 
However, neither the standard 3+1 equations nor the ADM equations are known to 
be well-posed, independently of the issue of constraint conservation. The ER and EB 
systems are well-posed, as are the equations describing conservation of the constraints 
in these systems. This subject is treated in more detail elsewhere [p^ . 
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